from charmicat import Cut
from cmix.parameters import paramdef

cuts = [Cut(name='pcm_pikpi_gt_2.7',
            desc='Select events with three-body momentum ' + \
                'greater than 2.7 GeV in the overall CM frame',
            pos_str='(pcm_pikpi > 2.7)',
            latex=r'$p^*(\pi_sK^+\pi^-) > 2.7$ GeV',
            ),
        Cut(name='m_pik_not_d0_28MeV',
            desc='Reject events where the K-pi+ invariant mass ' + \
                'is within 28 MeV of the D0 mass',
            pos_str='(abs(m_pik - D0_MASS) > 0.028)',
            condvars=paramdef,
            latex=r'$|M(K^-\pi^+) - M_{D^0}| > 28$ MeV',
            ),
        Cut(name='pos_kpi_llh_gt_0.5',
            desc='Identify K+ by atc_pid(K,pi) > 0.5',
            pos_str='(pos_kpi_llh > 0.5)',
            latex=r'$\mathcal{L}^+(K:\pi) > 0.5$',
            ),
        Cut(name='neg_kpi_llh_lt_0.9',
            desc='Identify pi- by atc_pid(K,pi) < 0.9',
            pos_str='(neg_kpi_llh < 0.9)',
            latex=r'$\mathcal{L}^-(K:\pi) < 0.9$',
            ),
        ]


